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Abstract 

A phenomenological model for the nucleon structure functions is presented. 
Visualising the nucleon as a cavity filled with parton gas in equilibrium and 
parametrizing the effects due to the finiteness of the nucleon volume, we 
obtain a good fit to the data on the structure function i*!*- The model then 
successfully predicts other unpolarized structure function data. 
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Recent experiments have revealed some remarkable features of the nucleon structure 
functions Ff'". Data on deep inelastic scattering of muons off proton and deuteron targets 

show that the quark sea in the nucleon is not flavor- symmetric, u{x) 7^ d{x); the Gottfried 
sum ^ Sg = I {Fi - F^) (dx/x), at = A GeV^, has the value 0.235 ± 0.026 compared 
to the usual quark model prediction of 1/3. This result has been confirmed by the observed 
asymmetry in Drell-Yan production of dileptons in pp and pn collisions Most notably, 
the HERA electron-proton scattering data 0] reveal a rapid rise of the proton structure 
function -Ff (2^) x decreases. Indeed over a wide range of small x, data from the various 
groups 0,3, for fixed Q^, are all well described by a single inverse power of x. Figure 1 is a 
log-log plot of the data on (the combination that enters So) versus x. We see that, 

for fixed Q^, the data fall on straight lines defined by 

^IM. = ±. , (0.0004 < X < 0.2). (1) 

For instance, at = 15 GeV^, the best-fit parameters are c = 0.229 ± 0.005 and m = 
1.22 ± 0.01.Q In Fig. 1, the straight-line extrapolations in the unexplored low-x region are 
our predictions based on Eq. (1). 

Global fits to the nucleon structure data involve parametrizing the various parton dis- 
tributions at some low and evolving them to higher relevant to observations. The 
fits so obtained have very high precision but contain several (typically ~15-20) arbitrary 
parameters and provide little physical insight into the structure of the nucleon. On the other 
hand, phenomenological models could give us some valuable clues into the physics of parton 
distributions in the nucleon. From this point of view the parton gas models or the statis- 
tical models of the structure functions have been quite interesting due to their intuitive 
appeal and simplicity. Bickerstaff and Londergan have provided a strong justification for 
the general philosophy of the statistical models and have also discussed limitations of other 



^For = 35, 120 GeV^ the power-law fits give, respectively, c = 
m = 1.26 ±0.01, 1.35 ±0.04. 
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0.229 ±0.006, 0.163 ±0.037 and 



models. 

We present here a phenomenological model for the unpolarized nucleon structure func- 
tions by using ideas from statistical mechanics. As a starting point, we assume that inside 
the nucleon, the valence quarks together with the sea quarks, antiquarks and gluons con- 
stitute a noninteracting gas in equilibrium. (It may be recalled that particles are treated 
as noninteracting even in the quark-parton model of Feynman.) This simple picture is then 
improved upon in two respects: (i) Finite-size (of the nucleon) corrections (FSC) to the 
statistical expression for the number of states per unit energy interval are taken into ac- 
count, and (ii) the resulting structure functions are evolved to the experimental values of 

by using the standard techniques in quantum chromodynamics (QCD). Each of these 
two effects is shown to play an essential role. To our knowledge, calculation of the structure 
functions taking into account FSC has never been reported in the literature. The model 
reproduces all unpolarized structure function data from x ~ 1 to x ~ 10~^ quite well. 

The Model 

We picture the nucleon (mass M) to consist of a gas of massless partons (quarks, an- 
tiquarks and gluons) in equilibrium at temperature T in a spherical volume V with radius 
R. We consider two frames, the proton rest frame and the infinite-momentum frame (IMF) 
moving with velocity — f (~ —1) along the common z axis. Our interest lies in the limit 
when the Lorentz factor 7 = (1 — t;^)^^/^ 00. The invariant parton number density in 
phase space is given by (quantities in the IMF are denoted by the index i) 

dn' ^ dn ^ gfjE) 
d^p^ d^r^ d^p d^r (27r)3 ' ^ ^ 

where g is the degeneracy {g = 16 for gluons and g = Q for g or g of a given flavor), {E, p) 
is the parton four-momentum and f{E) = {exp[f3{E — ^)] ± 1}^^ is the usual Fermi or Bose 
distribution function with (3 = T^^. In order to obtain the number distribution dn^ jdx in 
the Bjorken scaling variable x = pl/lMvy), we note that 

d^p' dV = 27rptrdpUMv-fdx)d^r/-f = 27iM^xdE dx d\ 
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and 

d^p d^r = AtxE^ dEd^r, 

for massless partons. For fixed x the parton energy E varies between the kinematic hmits 
{xM/2) and (M/2), where the lower hmit is attained when = and the upper hmit 
follows simply from energy-momentum considerations. (The kinematics will be described in 
detail in [l^]-) Hence drf/dx in the IMF is related to dn/dE and f{E) in the proton rest 
frame as follows 

dn' _ M^x /-A^/s dE dn _ gVM'^x r^l^ 

~ ^ J.M/2^dE ~ (27r)2 J,M/2 ^' 

The structure function ^2(2;) is given by 

The number distribution drf/dx in Eq. (3) vanishes linearly as x ^ (and also as x — 1) 
and leads to the behavior of the structure function F2{x) ~ at small x, which disagrees 
with the observations noted in Eq. (1). 

In an attempt to obtain the rise of F2{x) at small x, we now explore the effects aris- 
ing from the finiteness of the nucleon volume V. We note from Eq. (2) that dn/dE = 
gf{E)VE'^/ 271^, which is strictly valid in the large-volume limit, i.e., when the surface and 
curvature terms are negligible. We shall modify the model by incorporating these subleading 



terms. Various studies 0] of finite-size corrections (FSC) show that they depend on the 
particular equation of motion that is employed and are sensitive to the precise shape and 
size of the enclosure, the type of boundary conditions imposed on the wave function, and 
to the details such as whether the particles are strictly massless. Moreover, these studies 
invariably involve some simplifying assumptions and thus a blind use of their results cannot 
be justified in the present context. We feel more work needs to be done to fully understand 
the finite-size effects in a QCD bound state such as the nucleon. 

In keeping with the phenomenological nature of the model, we have incorporated the 



FSC by rewriting dn/dE in Eq. (3) as in [TT 



dn/dE = gf{E){VE'^/2TT^ + aR^E + hR), (4) 

and treating the numerical coefficients a and h in the surface and curvature terms as free 

parameters for reasons stated in the previous paragraph. 

The model described above is assumed to hold at a certain input momentum scale Qq; 

and if necessary can be evolved to higher by means of the standard techniques in QCD. 

To complete the statement of the model, we require the parton distributions to obey the 

following three constraints at the input scale. The constraints on the net quark numbers in 

the proton are n„ — = 2 and rid — n;i= 1, i-e., 

M2 /•! , fM/2 dE ( f dn\ f dn\ ] 

dx X -j^<\ — ]-{ — ]> = ria-ria. {a = u,d) (5) 



2 Jo JxM/2 ^2 t ydE ) ^ \dE^ 

Obviously, chemical potentials for heavy flavors are necessarily zero. As regards the third 
constraint, we assume that the longitudinal momentum fractions in the m, d flavors and the 
gluons add up to unity: 

^ rM/2 dE ( ( dn\ ( dn\ ( dn\ ( dn\ f dn\ ] 

-I "^^'Ue^ { (^)/ U)/ [^))-'- <«' 

The quark flavors s and c which are not introduced in Eq. (6) show up at higher as a 
result of QCD evolution. 

By interchanging the order of x and E integrations in Eqs. (5-6) and performing the 
x-integration analytically, we see that in order to keep the integrals finite, large powers of 
1/E are not allowed in the integrand. This means that the three terms displayed in Eq. (4) 
are the only ones that are allowed. Thus the model effectively has only two free parameters 
a and b. 

The temperature (T) and two chemical potentials /id) are not free parameters; we 
determine them by solving the three coupled nonlinear equations (5-6) by the Davidenko- 
Broyden method [0. The resulting values of T, fi^ and fid are such that the left and 
right hand sides of these equations agree with each other to typically one part in 10^. 
The parton distributions were evolved by means of the Dokshitzer-Gribov-Lipatov-Altarelli- 



Parisi equations [|13[ in leading order, taking the input scale Qq = M and Aqcd = 0.3 



GeV. Finally, the root-mean-square (rms) radius of the parton distribution was taken to be 



the same as the charge rms radius (p) of the proton; since p ~ 0.862 fm this yields 
R = ^sTs p = 1.11 fm. 

Results and Discussion 

Since the two arbitrary constants a and b in Eq. (4) are not known, we have determined 
them by fitting the deep inelastic scattering data on -Ffl^) = 15 GeV^ The 

results of our fit incorporating QCD evolution and FSC are shown by the solid curve in 
Fig. 2. Also shown for comparison in Fig. 2 are: (a) the (dot-dashed) curve labeled 'GAS' 
giving the prediction of the (unmodified) parton gas model which has no free parameters by 
virtue of the constraints, (b) the (dashed) curve labeled 'QCD' showing the effect of QCD 
evolution on the gas model, and (c) the (dotted) curve labeled 'FSC showing a fit to the 
data when only the FSC are introduced in the gas model. For the solid curve in Fig. 2, the 
fitted values of the two parameters are a = —0.400 and b = 0.475, and the corresponding 
temperature and chemical potentials are T = 72 MeV, p„ = 162 MeV and p^ = 81 MeV.Q 

As a test of the model, we show in Fig. 3 the prediction (solid curve) for the difference 
[F2 (x) — F2{x)]. Also shown for comparison is the result (dashed curve) based on the 
parametrization of Gliick et al. The agreement with the NMC data is reasonable. 

As for the Gottfried sum 5*^, we have 

1 2 M2 /-i r^^/^ dE ( ( dn\ ( dn\ 



The inequality 6*0 < | is thus a result of having in the proton, more valence u quarks than 
valence d quarks, (ra^ — «„) > {rid — n^), implying that p^ > p^ (see Eq. (5)) and hence the 
integral in So is positive. Our model predicts at = 4 GeV^, the value Sg = 0.22 which 



^It is amusing to note that the values of a and b determined by us are close to the values o = — 1/2 
and b = 3/(27r) which follow from one of the expressions for dn/dE given by Morse and Ingard 



|rT|; substitute c = 1, 1/ = E/{2t:),A = AttR'^ and = £y = £^ = 2R in dNob in their Eq. (9.5.12). 
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is consistent with the experimental value So = 0.235 ± 0.026. The natural explanation of 
the violation of the Gottfried sum rule is an attractive feature of our general framework. 

In addition, the model is in excellent agreement with data on the gluon distribution g{x). 
It also reproduces well the ratio F2{x)/F2{x), the ratio u{x)/d{x) at < x >= 0.18 which has 
been deduced to be about 0.51 by the NA51 collaboration [Q, the longitudinal momentum 
fraction carried by the charged partons, etc. These and other predictions of the model, on 
the quark and antiquark distributions q{x), q{x), q.o{x) = q{x) — q{x) for various flavors, etc. 
will be discussed in detail elsewhere |llO[| . 

Now we briefly describe the salient features of some of the recent calculations of the 
nucleon structure functions, which use ideas from statistical mechanics. It has generally 
been difficult so far to model F2{x,Q'^) such that it is nonzero only in x = [0, 1], and is 
consistent with the two number constraints, the momentum constraint, the Gottfried sum 
and the observed low-x rise. Mac and Ugaz [7a] calculated first-order QCD corrections 
to the statistical distributions and obtained a crude but reasonable agreement with -^^(x) 
data for x ^ 0.2. The momentum constraint was not imposed and the fitted value of the 
proton radius was 2.6 fm. Cleymans et al. [7b] used the framework of the finite temperature 
quantum field theory. They considered 0{as) corrections to the statistical distributions 
and obtained a good fit to the F2{x) data for x > 0.25. They also calculated the ratio 
c"l / <^T in this region; it was a factor of 6 above the experimental value. Bourrely et al. [7c] 
considered polarized as well as unpolarized structure functions and presented a statistical 
parametrization (with eight parameters) of parton distributions in the IMF. Their framework 
allowed chemical potential for quarks as well as for gluons. The number constraints were 
not satisfied very accurately. QCD effects were not considered. Parton distributions were 
nonzero for x > 1. xq{x) vanished as x — > and so it was not possible to reproduce the fast 
increase of the antiquark distributions for x < 0.1. Bourrely and Soffer's [7d] approach was 
similar to that in [7c] . By incorporating QCD evolution of parton distributions and allowing 
the antiquark chemical potential to depend on x, they were able to reproduce the HERA 
data on . Recently Buccella et al. [7e] have introduced in this model a so-called liquid 
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term to reproduce the low-x behaviour of structure functions. 

In conclusion, it is noteworthy that the apphcation of ideas of statistical mechanics to the 
point constituents of the nucleon can provide a simple description of all the observed features 
of the unpolarized nucleon structure functions down to the lowest x values so far explored. 
The model presented here has two free parameters which arise from our treatment of the 
finite-size corrections. It is hoped that the success of the model would provide a stimulus to 
further studies of the finite-size effects in a QCD bound state such as the nucleon. 
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FIGURES 

FIG. 1. Log- log plot of the proton structure function data. Experimental data are from Refs. 

[4,5]; the error bars show statistical and systematic errors combined in quadrature. The straight 
lines are our fits described in Eq. (1), and are labeled by = 15, 35, and 120 GeV^. Numbers 
have been scaled by the factors shown in parentheses for convenience in plotting. 

FIG. 2. Proton structure function (^) — GeV^. Data points are as in Fig. 1. Solid 

curve is our best fit to the data and includes both FSC and QCD. Also shown for comparison are: 
the (dot-dashed) curve labeled 'GAS' giving the (unmodified) gas model prediction, the (dashed) 
curve labeled 'QCD' showing the QCD-evolved gas model, and the (dotted) curve labeled 'FSC' 
which is a fit to the data when finite-size corrections are included in the gas model (without QCD). 

FIG. 3. Difference (-Fg ~ ^2) versus x, at = 4 GeV^. Experimental data arc from Rcf. [1]; 
errors are statistical only. Solid curve is the prediction of our model. Dashed curve is based on the 
parametrization of Gliick et al. [6] . 
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